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Abstract 

Avalanche  gains  >80  are  observed  in  InAsyPl-y  (0.1<y<0.3)  grown  on  p-type  InP  substrates  with 
lattice  strains  <1%.  The  devices  have  a  primary  dark  current  >200  pA.  The  InAsyPl-y  layers  have  the 
same  lattice  parameter  as  compositions  of  InxGai.xAs  which  absorb  wavelengths  <2.1|im,  for  use  in 
avalanche  photodiode  imaging  arrays  operating  in  the  mid-IR  spectral  range. 
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Avalanche  Gain  in  InASyPi.y  Photodetectors  for  Mid-IR  Focal  Plane  Arrays 


D.  S.  Kim,  S.  R.  Forrest,  G.  H.  01sen@,  M.  J.  Lange@,  R.  U.  Martinelli#  and  N.  J.  Di  Giuseppe# 

Advanced  Technology  Center  for  Photonics  and  Optoelectronic  Materials 

Princeton  University 
Princeton,  NJ  08544 

Very  high  performance,  integrated  focal  plane  arrays  for  use  in  imaging  applications  have  been 
demonstrated  using  lattice-matched  Ino.53Gao.47As/InP  p-i-n  photodiodes  as  the  detection  element. 

Unfortunately,  these  devices  are  sensitive  only  to  wavelengths,  X<1.7pm,  limiting  their  utility  for 
many  moisture  sensing,  LIDAR  and  thermal  measurement  applications.  More  recently,  we  have 
demonstrated  integrated  "three-color"  detectors  employing  InxGai_xAs  (x  >  0.53)  photodiode  arrays 

which  have  optical  sensitivity  to  X  =  2.6pm  [1].  For  very  low  signal  levels,  it  is  useful,  however,  to 
extend  the  sensitivity  obtained  using  a  simple  p-i-n  detector.  One  possibility  for  increasing  long 
wavelength  focal  plane  sensitivity  is  to  use  avalanche  photodiodes  (APDs).  In  this  work,  we  describe 
initial  success  in  the  effort  to  develop  such  a  device.  We  report  the  observation  of  large  avalanche 
gains  (40-80)  in  the  presence  of  small  primary  dark  currents  (-500  pA  at  90%  of  the  breakdown 
voltage,  Vg)  in  InAsyPl-y,  where  0.1  <  y  <  0.3  was  chosen  to  make  this  material  lattice-matched  to 
narrow  band  gap  (~0.5eV)  InxGai-xAs  (0.6<x<0.7).  These  results  represent  the  first  significant  step 
in  realizing  APDs  for  integrated  focal  plane  arrays  useful  in  the  mid-IR  wavelength  region. 

The  device  structure  to  be  used  for  mid-IR  detection  is  shown  schematically  in  Fig.  la.  The 
structure  is  a  separate  absorption  and  multiplication  region  APD  (or  SAM-APD)  which  has  proven 
effective  in  eliminating  tunneling  leakage  current  that  occurs  in  narrow  band  gap  semiconductors  such 
as  In0.53Gao.47As  when  in  the  presence  of  the  high  applied  electric  fields  required  for  avalanche 

multiplication.  To  make  such  a  device  sensitive  to  X  >1.65pm,  additional  In  must  be  added  to  InxGai- 
xAs.  Since  this  results  in  an  increase  in  the  lattice  constant,  this  material  is  no  longer  lattice-matched  to 
the  large  band  gap  InP.  Lattice-mismatch  generates  misfit  dislocations  at  the  interface  between  the 
contacting  layers,  which  act  as  recombination  centers  that  increase  dark  current,  and  serve  as 
nucieation  sites  for  microplasmas.  To  circumvent  this  problem,  As  is  added  to  the  InP  multiplication 
region,,  such  that  InAsyPj.y  serves  as  the  large  band  gap  material.  The  As  concentration  is  chosen 
such  that  this  layer  will  ultimately  have  the  same  lattice  constant  as  the  absorbing  InxGa].xAs.  In  this 
way,  dislocations  are  confined  to  the  mismatched  InP/InAsyPi-y  interface,  while  being  minimized  at 
the  lattice-matched  InAsyPi-y/InxGai.xAs  heterojunction. 

To  determine  the  gain  characteristics  of  InAsyPi-y,  the  structure  in  Fig.  lb  was  fabricated. 
Undoped,  3pm  thick  n"-InAsyPi-y  was  grown  by  MOCVD  directly  onto  a  p-type,  Zn:InP  substrate. 
During  growth  at  700°C,  Zn  from  the  highly  doped  (>5x10^  cm_3)  substrate  diffuses  into  the 
InAsyPi.y  layer,  forming  the  p-n  junction.  This  places  the  p-n  junction  (which  is  the  location  of  the 
highest  electric  fields)  beyond  the  extent  of  the  misfit  dislocations  that  form  at  the  InAsyPi_y/InP 
heterojunction.  Similar  results  are  obtained  for  devices  grown  by  hydride  vapor  phase  epitaxy. 
Following  growth,  a  50  pm  diameter  Au:Ge:Au  n-contact  was  deposited  onto  the  surface  of  an 
approximately  150  pm  diameter  mesa  structure,  formed  by  etching  in  1%  Br:MeOH.  Contact  alloying 
was  achieved  by  rapid  thermal  annealing. 

Three  wafers  were  grown  and  tested  for  this  study,  with  their  parameters  listed  in  Table  1.  In 
the  far  right-hand  column  of  the  Table  is  shown  the  cut-off  wavelength  of  the  composition  of  InxGai„ 
xAs  lattice-matched  to  InyAsj.yP  in  the  adjacent  column.  To  measure  gain,  the  devices  were 
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illuminated  from  the  top  surface  using  633  nm  light  from  a  He:Ne  laser,  or  from  the  bottom  surface 
using  1064nm  light  from  a  Nd-YAG  laser.  The  InP  substrate  is  transparent  to  the  1064nm  radiation. 
Light  was  introduced  using  a  single  mode  fiber  scanned  from  edge  to  edge  of  the  mesa  to  ensure  that 
the  gain  of  the  devices  was  uniform.  Since  the  633nm  light  is  strongly  adsorbed  at  the  InAsyPi-y 
surface,  a  measurement  of  gain  due  to  top  illumination  yields  the  hole  multiplication  factor,  Mp. 
Similarly,  illumination  through  the  substrate  side  provides  the  electron  multiplication  factor,  Mn. 

The  dark-current  is  relatively  voltage  independent  at  V<  25V,  increasing  monotonically  toward 
breakdown  at  voltages  greater  than  40V.  Primary  dark  currents  as  low  as  200  pA  for  wafer  A  are 
observed  at  0.5Vb  =  25V,  whereas  the  more  highly  strained,  lower  band  gap  wafer  C  has  a  dark 
current  of  100  nA  at  this  same  voltage.  In  the  case  of  wafer  A,  the  dark  current  increases  to  only  InA 
at  0.9Vb.  These  results  suggest  that  the  dark  current  is  increasing  primarily  as  a  result  of  avalanche 
gain  rather  than  the  result  of  parasitic  effects  such  as  edge  breakdown  or  tunneling.  The  photocurrent 
and  multiplication  characteristics  of  the  wafer  C  are  shown  in  Fig  2.  The  photocurrent  remains 
constant  until  approximately  0.6Vb,  at  which  point  M  increases  toward  its  maximum  value  at  Vb- 
From  these  data,  we  observe  that  gains  as  high  as  70  -  80  are  obtained  for  wafers  B  and  C,  and  from 

30  to  40  for  wafer  A.  Similarly,  at  all  voltages  and  for  all  wafers,  we  find  that  Mp  =  Mn.  Variations  in 
both  the  maximum  values  of  M,  and  between  Mn  and  Mp  from  wafer  to  wafer  are  not  significant. 

The  electric  fields  at  breakdown  are  measured  to  be  1-L5xl05  V/cm,  or  roughly  four  times 
smaller  than  is  required  in  conventional,  III-V  based  APD  materials  such  as  InP.  The  probable  origin 
of  the  secondary  carriers  is  ionization  of  the  deep  levels  arising  from  the  misfit  defects,  which  would 
also  account  for  the  nearly  equal  values  obtained  for  both  Mn  and  Mp  in  these  diodes.  To  test  this 
conclusion,  we  also  measured  the  gain  as  a  function  of  light  intensity  for  633nm  wavelength  light 
incident  via  the  top  surface  (Fig.  3).  The  gain  is  a  strong  function  of  light  intensity,  with  Mp  at  46V 
decreasing  from  a  maximum  value  of  40  at  an  input  power  of  350  nW,  to  only  9  at  6  (iW.  This  is 
consistent  with  the  hypothesis  that  gain  occurs  via  the  ionization  of  a  limited  number  of  defect  states. 
Hence,  the  gain  decreases  at  relatively  high  input  power  levels  where  a  large  fraction  of  the  defects  are 
ionized  at  any  given  time. 

Preliminary  results  for  a  full  SAM-APD  structure  operating  in  the  mid-IR  range  will  be 
presented. 

In  conclusion,  we  have  demonstrated  that  low  primary  dark  currents  and  high  avalanche  gains 
can  be  achieved  in  lattice-mismatched  InAsyPi-y  (0.1<y<0.3)  grown  without  compositional  grading  on 
p-type  InP  substrates,  where  the  InAsyPi-y  strain  is  as  high  as  1%.  This  range  of  InAsyPj.y  material 
is  lattice-matched  to  InxGai-xAs  (0.6<x<0.7),  which  can  absorb  light  at  wavelengths  as  long  as  2.1 
(im.  Hence,  these  results  strongly  suggest  that  InAsyPi-y  can  be  used  as  the  multiplication  layer  in 
SAM-APD  imaging  arrays  useful  for  sensing  applications  in  the  mid-IR.  We  have  also  found  that  the 
ionization  coefficients  of  electrons  and  holes  are  approximately  equal  in  the  compositions  of  InAsyPj. 
y,  and  over  the  range  of  electric  fields  investigated. 

The  authors  acknowledge  the  support  of  Ballistic  Missile  Defense  Organization  and  the  Office 
of  Naval  Research  (A.  M.  Goodman  and  Y.  S.  Park). 

Table  1:  Parameters  of  InAsyPi-y  Wafers 


Wafer  ID 

y 

lattice  parameter 

(A) 

strain  (%) 

A-gJ  InASyPj-y 
(Hm) 

A,g;  InGaxAsi-x 
(tim) 

A 

0.12 

5.8947 

-0.43 

1.02 

1.87 

B 

0.25 

5.9199 

-0.71 

1.15 

2.04 

C 

0.29 

5.9267 

-0.97 

1.20 

2.11 

[1]  D.  S.  Kim,  S.  R.  Forrest,  M.J.  Lange,  G.  H.  Olsen  and  M.  J.  Cohen,  Photonic  Tech.  Lett.  f>,  235 
(1994). 
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An  8  Channel  Digitally  Tunable  Transmitter 


C.H.  Joyner,  M.  Zirngibl  and  J.C.  Centanni 
AT&T  Bell  Laboratories 
Crawford  Hill  Laboratory 
Holmdel,  NJ  07733 


Abstract: 

We  demonstrate  an  8  frequency  (1.567  to  1.580  urn  )  digitally  tunable  laser  based  on 
a  waveguide  grating  router  with  an  electroabsorption  modulated  (-10dBA/)  output.  A 
unique  output  coupler  design  delivers  an  ASE  noise  suppression  of  -55  dB/A. 
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An  8  Channel  Digitally  Tunable  Transmitter  with  Electroabsorption 
Modulated  Output  by  Selective  Area  Epitaxy 

C.H.  Joyner,  M.  Zirngibl  and  J.C.  Centanni 
AT&T  Bell  Laboratories 
Crawford  Hill  Laboratory 
Holmdel,  NJ  07733 


I.  INTRODUCTION 

Wavelength-division-multiplexed  (WDM)  lightwave  communication  systems  are 
receiving  increasing  attention  for  long  distance  as  well  as  local  applications.  Many 
such  multi-frequency  lasers  have  been  proposed  [1-4],  but  few  have  been  integrated 
with  modulators  to  become  transmitters.  In  this  paper  we  present  an  intregrated  WDM 
transmitter  transmitter  with  precise  channel  spacings  in  which  any  of  8  wavelengths 
may  be  selected  and  modulated  at  a  single  output  port. 

II.  DESIGN  AND  FABRICATION 

A  schematic  diagram  of  the  transmitter  is  shown  in  figure  1.  The  laser  cavity  is  similar 
to  our  previously  reported  design  [5].  The  novel  feature  of  the  current  design  arises 
from  the  realization  that  light  emerging  from  the  waveguide  grating  arms  radiates  into 
first  order  and  second  order  grating  modes.  The  primary  mode  is  always  aligned  with 
the  amplifier  completing  the  left  side  of  the  laser  cavity,  while  the  secondary  mode 
represented  an  unwanted  loss  in  past  designs.  In  the  new  design,  the  free  space 
region  now  includes  an  output  port  that  captures  second  order  diffracted  light  from  the 
waveguide  grating  ( ref.  M.  Zirngibl,  "Wavelength  grating  router  with  output  couplers", 
U.S.  patent  application).  This  port  serves  as  an  output  coupler  from  the  main  lasing 
cavity  without  increasing  the  loss  of  the  lasing  cavity.  The  output  port  is  coupled  to  an 
on-chip  electroadsorption  modulator  that  is  terminated  on  an  AR  coated  facet. 

The  cross-section  of  our  device  is  shown  in  figure  2.  It  is  similar  to  our  previously 
reported  structure  [6],  except  that  the  laser  active  region  and  modulator  are  grown  by 
selective  area  epitaxy  (SAE). 

III.  RESULTS  and  DISCUSSION 

The  device  was  tested  under  DC  conditions  in  chip  form.  With  cleaved  facets,  the 
laser  had  a  cw  threshold  of  60-70  mA  with  this  current  split  equally  between  the 
amplifiers  at  either  end  of  the  cavity.  The  output  power  measured  directly  at  the  left 
cleaved  laser  facet  was  -10  dBm  into  a  lensed  single  mode  fiber,  with  an  applied 
current  of  100  mA.  The  power  from  the  modulator  port  was  -16  dBm  under  the  same 
conditions.  The  modulator  gave  an  extinction  ratio  of  -10  dBA/  at  a  bias  of  -5V. 

Passive  waveguide  tests  with  the  amplifiers  removed  indicate  that  the  coupling  to  the 
modulator  port  is  approximately  -4  dB.  The  upper  half  of  Figure  3  shows  spectra  in 
which  each  of  the  channels  was  powered  separately,  monitored  at  the  modulator 
facet,  and  superimposed.  The  relative  channel  spacing  is  extremely  precise.  Note, 
however,  that  channels  1  and  8  also  have  significant  power  in  the  next  free  spectral 
range  of  the  router.  The  lower  half  of  figure  3  shows  the  output  of  channel  5  observed 
through  the  modulator  port  on  a  more  sensitive  power  scale  with  lA  bandwidth 
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resolution.  Note  that  the  amplified  spontaneous  emission  (ASE)  noise  is  suppressed 
by  -55dB  when  measured  with  1  A  resolution,  which  shows  another  advantage  of  our 
output  coupler  design.  The  signal  is  filtered  by  one  additional  pass  through  the  router 
before  it  reaches  the  modulator.  If  the  output  fiber  were  coupled  to  the  amplifier  port, 
noise  suppression  would  typically  be  only  -40  dB  /  A.  The  contributions  from  other 
pass  bands  of  the  router,  seen  left  and  right  of  the  main  peak,  are  generally  down  by 
-35dB  for  channels  2  through  7. 

The  reason  for  the  anomalous  behavior  of  channels  1  and  8  is  the  nature  of  the  SAE 
mask  used  to  create  the  band  gap  shifts.  SAE  is  sensitive  to  the  filling  factor  (ratio  of 
area  masked  to  that  unmasked)  in  the  immediate  vicinity  of  the  unmasked  region 
undergoing  epitaxial  growth.  The  amplifiers  for  channels  1  and  8  were  at  either  outer 
edge  of  the  amplifier  array.  Thus,  seeing  neighbors  on  only  one  side,  they 
experienced  less  of  an  increase  in  quantum  well  thickness  and  therefore  a  shorter 
wavelength  bandgap.  The  main  determining  factor  in  the  selection  of  which  bandpass 
the  cavity  will  lase  upon  is  the  peak  of  the  ASE  gain  spectrum.  This  difficulty  can  be 
easily  overcome  by  using  dummy  amplifier  sections  flanking  the  outermost  ones. 

IV.  CONCLUSIONS 

We  have  demonstrated  an  8  frequency  digitally  tunable  laser  based  on  a  waveguide 
grating  router  with  an  electroabsorption  modulated  output.  The  channel  wavelengths 
range  from  1.567  to  1.580  pm,  useful  for  WDM  communications  networks.  The 
channel  spacings  of  (1.6nm)  200  GHz  are  extremely  precise.  The  DC  modulator 
extinction  ratio  is  -10  dBA/  with  a  -5V  bias.  A  unique  output  coupler  design  delivers  an 
ASE  noise  suppression  of  -55  dB/A  . 
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Abstract 

We  report  the  first  polarization-independent  monolithic  WDM  detector.  8  x  2nm  channels, 
from  1.544^m  to  1.558/rm  wavelength,  were  demultiplexed  and  detected  by  the  compact  3.7mm  x 
4.0mm  device  without  measurable  polarization-sensitivity. 

Summary 

Introduction: 

Multiple-wavelength  detection  is  needed  for  emerging  wavelength  division  multiplexed 
(WDM)  fiber  networks.  One  approach  being  researched  is  the  integration  of  the  wavelength 
demultiplexing  and  detection  functions  onto  a  single  chip  of  InP-based  material.  This  avoids  the 
problem  of  connecting  an  array  of  fibers  to  detectors  and  allows  the  possibility  of  monolithic 
integration  of  array  amplifiers  to  form  a  single-chip  WDM  receiver.  High-speed" detection  using 
small  area  waveguide  detectors  is  also  readily  achievable  [1], 

One  of  the  most  significant  impediments  to  this  approach,  however,  is  that  InP-based 
wavelength  demultiplexers  are  generally  birefringent,  dispersing  TE  and  TM  signals  differently. 
The  InP-based  WDM  detectors  reported  to  date,  based  on  either  planar  reflection  grating 
demultiplexers  [2,  3]  or  arrayed  waveguide  demultiplexers  [4,5],  have  had  TE  -  TM  signal 
separations  ranging  from  0.6  to  2.6  times  the  channel  spacing.  This  is  clearly  unacceptable  for 
most  fiber  systems,  where  the  polarization  state  of  the  input  light  is  undetermined. 

In  this  paper  we  report,  for  the  first  time,  a  polarization-independent  InP-based  WDM 
detector.  It  is  based  on  an  arrayed  waveguide  demultiplexer  that  employs  raised-stripe  waveguides 
designed  to  be  free  of  modal  birefringence.  We  have  integrated  this  with  an  array  of  waveguide  p- 
i-n  diodes  and  demonstrate  polarization-independent  detection  for  8  channels  separated  by  2nm. 

Design  and  fabrication: 

The  array  demultiplexer  employed  InGaAsP(Eg=1.0^m)  raised-stripe  guides  on  InP. 
Calculation  indicated  that  1.9/rm  thickness  and  3/<m  stripe-widths  give  n(TE)=n(TM).  A  compact 
3.7mm  x  2.6mm  demultipexer  was  designed  integrated  with  detectors  employing  the  “hybrid” 
vertical  /  butt-coupled  scheme  [ l ].  A  diagram  of  the  guide  and  detectors  structures  is  given  in 
Figure  1.  Detector  layout  is  shown  in  Figure  2.  The  die  size  of  the  detector  was  3.7mm  x  4.0mm. 

The  layer  structure  was  grown  by  LP-OMCVD.  First,  the  detector  layers  (1.0/rrn  InGaAs 
absorption  layer;  0.6/rm  p-1.3Q/p-InGaAs  contact  layer)  were  grown  on  n-InP,  above  an  0.6/<m- 
thick  n-doped  portion  of  the  1.0Q  waveguide  layer.  Following  mesa-definition  of  the  detectors,  a 
further  1.3/<m  of  s.i.-l.OQ  was  grown  to  complete  the  waveguiding  layer  and  provide  isolation. 
Waveguide  stripes  were  then  etched  by  H2/CH4  R.I.E.,  ~0.2;<m  into  the" underlying  InP.  Detector 
contacts  were  made  by  Ti/Pt/Au  mesa-roof  metallisation  and  Ni/Ge/Au/Ni/Au  back-wafer 
metallisation.  The  wafer  was  then  thinned  and  cleaved  to  yield  single  WDM  detector  chips.  No 
A.R.  coating  was  applied  to  the  input  facet. 

Array  demultiplexers  both  with  integrated  detectors  and  without  detectors  were  fabricated. 
A  tunable  laser  and  polarization  controller  was  used  with  a  lensed  fiber  to  provide  TE  orTM  input 
and  the  detector  responses  were  measured  on-chip  at  zero  bias. 

Results  and  discussion: 

The  response  of  one  detector  for  incident  TE  and  TM  signals  is  shown  in  Figure  3.  Within 
the  measurement  accuracy,  the  detector  responses  for  TE  and  TM  polarized  light  were  identical;  the 
spectral  variation  forTE  and  TM  light  and  the  absolute  channel  responses  are  indistinguishable. 


Detector  response 
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The  small  difference  in  signal  levels  for  the  different  channels  (~1.4dB  in  worst  case)  is  ascribed  to 
defects  in  the  feed-guides  leading  to  the  detectors. 

The  wavelength  channel  spacings  were  found  to  be  exactly  2nm  (±0.06nm  r.m.s.,  on  a 
spectral  step  of  O.lnm),  with  no  noticeable  deviation  across  the  8-channel  span;  the  8-channel- 
span  measured  14.0nm,  within  the  O.lnm  accuracy  of  the  scans. 

The  detectors  recorded  a  passband  FWHM  of  ~0.7nm,  in  line  with  expectations,  and 
identified  a  passband  /  channel  spacing  ratio  of  35%.  Cross-talk  between  channels  was  typically  - 
16dB,  though  somewhat  higher  for  adjacent  channels,  at  ~-13dB.  Similar  values  were  observed 
for  demultiplexers  that  did  not  have  photodetectors  integrated,  showing  it  to  be  optical  in  nature. 
The  background  value  is  higher  than  the  best  reported  figures  for  raised-stripe  array  demultiplexers 
(<-25dB,  [6]),  and  is  possibly  associated  with  waveguide  phase  fluctuations  arising  from  guide- 
width  variations.  This  is  consistent  with  typical  demux  excess  losses  of  6.5-7.5dB,  which  is  also 
higher  than  the  value  of  ^4-5dB  reported  in  that  reference.  (Some  loss  will  also  be  associated  with 
the  planar/stripe-guide  junctions.)  Total  on-chip  loss  of  the  detector  shown  was  9-10dB,  with  an 
in-coupling  loss  of  3.9dB  (no  AR-coating).  Detector  external  responsivities  were  0.05-0.065A/W. 

Conclusions: 

We  have  demonstrated,  for  the  first  time,  an  integrated  polarization-insensitive  WDM 
detector.  It  comprised  an  arrayed  waveguide  demultiplexer  employing  raised-stripe  1.0Q  InGaAsP 
guides  and  a  p-i-n  waveguide  detector  array.  8-channel  WDM  detection  was  demonstrated  at 
precisely  2nm-spaced  channels  without  measurable  polarization  sensitivity.  On-chip  losses  were 
9.0-10.0dB  and  channel  cross-talk  was  -13dB  to  -16dB. 
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Fig.  3:  TE  response  (dashed)  and  TM  response  (solid). 
Also, detail  for  one  channel. 
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ABSTRACT 

A  four  channel  Multi-Grating  Filter  (at  1.55  microns)  with  polarization-independent 
channel  positions  and  two  integrated  amplifiers  has  been  fabricated  in  InGaAsP/InP. 
A  packaged  device  has  no  fiber-to-fiber  loss  and  worst  case  crosstalk  below  -18  dB. 
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One  of  the  key  component  for  the  practical  real¬ 
ization  of  Wavelength  Division  Multiplexing  (WDM) 
transmission  systems  is  a  tunable  optical  filter  or 
demultiplexer  that  allows  the  selection  of  the  de¬ 
sired  channel.  Reported  devices  suitable  for  integra¬ 
tion  include  filters  based  on  active  and  passive  DBR 
or  DFB  structures,  (multiple)  Mach-Zender  interfer¬ 
ometers  and  co-directional  grating-assisted  couplers 
(see  references  in  [1]).  Promising  approaches  include 
acousto-optic  filters  [2]  and  demultiplexers  based  on 
waveguide  arrays  [3,  4]  or  etched  diffraction  gratings 
[5,  6,  7].  We  have  taken  yet  a  different  approach  with 
a  device  called  a  Multi-Grating  Filter  (MGF)  [1,  8,  9]. 
The  MGF  acts  as  a  demultiplexer  by  allowing  only 
one  channel  (or  eventually  several  channels)  of  a  dis¬ 
crete  set  to  be  transmitted  and  blocking  the  other 
ones,  but  it  acts  as  a  filter  in  that  the  selection  of  the 
transmitted  channel  is  controlled  electrically. 

To  be  useful  in  a  real  system,  these  filters  and  de- 
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Figure  2  : Normalized  transmission  spectrum  of  a 
packaged  MGF  (# 7 )  using  the  (unpolarized)  spon¬ 
taneous  emission  from  an  Erbium  doped  Fiber  Am¬ 
plifier  (amplifier  currents  were  35  mA  and  15  mA). 
[Add  1.6  nm  to  get  the  correct  wavelengths] 
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Figure  1:  Structure  of  the  four-channel  polarization 
independent  MGF  with  two  integrated  amplifier  sec¬ 
tions. 

multiplexers  have  to  be  polarization  insensitive.  The 
MGF  devices  reported  earlier  [1,  8,  9]  had  about 
4.5  nm  difference  between  the  channel  wavelengths 
for  the  TE  and  the  TM  mode.  In  this  paper,  we 
present  a  Multi-Grating  Filter  with  a  new  waveguide 
design  which  eliminates  the  TE/TM  offset.  Pack¬ 
aged  devices  with  no  fiber-to-fiber  loss  and  worst  case 
crosstalk  below  -18  dB  have  been  fabricated. 

As  described  previously  in  [1,  8,  9],  the  MGF  uses 
several  Bragg  gratings  in  series  in  a  waveguide  (see 
Figure  1).  Each  grating  center  wavelength  corre¬ 
sponds  to  one  channel  wavelength,  so  that  all  chan¬ 
nels  are  blocked.  A  channel  is  selected  by  injecting 
current  in  the  corresponding  grating,  which  changes 
the  center  wavelength  and  allows  the  transmission  of 
that  channel.  The  center  wavelength  of  a  grating  Ac  is 
given  by  \c~2ncj j  where  A  is  the  (physical)  grat¬ 
ing  period  and  nejj  the  effective  index.  The  different 
gratings  are  fabricated  with  different  A  and  the  cur¬ 
rent  injection  changes  nejf  through  the  plasma  effect 
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Figure  3:  Packaged  MGF  #7  :  transmission  of  channel  f  (J552Jf  nm)  as  a  function  of  the  LA  current  and 
input  power  in  the  fiber.  SA  is  biased  at  transparency  for  that  channel.  MIN  and  MAX  values  for  varying 
input  polarizations  are  plotted. 


of  the  injected  carriers  (see  [10]).  The  TE  and  TM 
modes  will  have  the  same  center  wavelengths  if  they 
have  the  same  nej / . 

This  leads  to  a  waveguide  design  with  a  nearly 
square  cross-section.  But  there  are  several  other  con¬ 
straints:  single  mode  waveguide,  large  grating  cou¬ 
pling  coefficient  k  (around  100  cm"1),  sufficient  cen¬ 
ter  wavelength  shift  with  current  injection  (2  to  3 
nm),  waveguide  as  wide  as  possible  (>  0.5  /zm)  to 
allow  easy  fabrication,  reasonable  tolerances.  Using 
a  vectorial  finite  element  program  to  compute  the 
effective  indices,  we  arrived  at  the  design  shown  in 
Figure  1.  Using  the  refractive  indices  from  [10]  (at 
300K,  1.555  /im),  calculations  give  ntjj—  3.223,  k~ 
107  cm"1,  TE-TM  difference^  0.05  nm.  The  most 
difficult  dimension  to  control  is  the  waveguide  width. 
Calculations  show  that  going  from  0.7  /zm  to  0.8  /zm 
increases  the  TE-TM  offset  to  0.25  nm  and  k  to  117 
cm"1,  while  the  center  wavelength  shifts  by  3.2  nm. 
Thus,  the  waveguide  width  should  be  controlled  to 
better  than  0.05  /im,  which  is  possible  with  electron- 
beam  lithography  and  dry  etching. 

The  devices  (shown  in  Figure  1)  are  fabricated  on 
a  n-type  InP  substrate  using  low  pressure  MOVPE 
(Metal  Organic  Vapour  Phase  Epitaxy)  for  all  the 
epitaxial  growths.  First  the  layers  for  the  amplifier 
sections  are  grown:  lower  n-InP  cladding,  0.6  /zm  n- 
InGaAsP  (1.1  /zm)  waveguide  layer,  50  nm  n-InP,  50 
nm  n-InGaAsP  (1.4  /zm)  grating  layer,  0.1  /z m  n-InP, 


0.2  /zm  undoped  InGaAsP  (1.58  /zm)  active  layer  and 
top  p-InP  cladding.  All  n  dopings  are  1018  cm"3 
and  the  p  doping  is  6  x  1017  cm-3.  Then  the  grat¬ 
ing  regions  are  etched  down  to  the  grating  layer  and 
the  gratings  are  formed  by  electron-beam  lithogra¬ 
phy  and  etching  in  the  grating  layer.  Then,  50  nm 
n-InP  is  regrown  over  the  gratings,  followed  by  a  0.5 
/zm  undoped  InGaAsP  (1.2  /zm)  top  waveguide  layer 
and  a  top  p-InP  cladding  layer.  The  0.7  /zm  wide 
waveguide  is  then  defined  by  electron-beam  lithog¬ 
raphy  and  etched  using  reactive  ion  etching.  After 
etching  of  the  top  p  cladding  between  sections  (for 
electrical  isolation),  semi-insulating  (Fe  doped)  InP 
is  selectively  grown  for  current  confinement.  Finally, 
the  top  p  contact  layers  are  grown  and  metalization  is 
done.  Then  the  contact  layers  are  removed  between 
the  sections.  Each  grating  section  is  400  /zm  long  and 
the  amplifiers  are  450  /zm  (LA)  and  250  /zm  long  (SA) 
respectively.  The  waveguides  terminate  20  /zm  from 
the  anti-reflection  coated  facets  (to  reduce  reflection). 
The  filter  was  packaged  the  same  way  as  in  [9],  with 
a  Peltier  element,  thermistor  and  two  lens  couplers, 
except  that  we  used  normal  fibers  and  no  isolator. 

Using  polarization  maintaining  fiber  tapers  to  cou¬ 
ple  light  in  and  out  of  the  device  and  the  spontaneous 
emission  from  an  EDFA  as  the  light  source,  we  mea¬ 
sured  less  than  0.2  nm  difference  in  center  wavelength 
between  TE  and  TM  (before  packaging).  Some  de¬ 
vices  had  no  measurable  difference  (<  0.1  nm).  An 
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Open 

Wavelength 

Gain 

[dB] 

On-off 

Crosstalk  [dB 

from  channel 

[  Transparency 

channel 

[nm] 

max 

min 

ratio  [dB] 

1 

2 

3 

4 

SA  [mA] 

LA  [mA] 

i 

1542.1 

+1.7 

-0.7 

17.0 

- 

-22.4 

-22.7 

-25.0 

5.2 

9.8 

2 

1545.8 

+2.8 

-1.3 

18.7 

-18.5 

- 

-20.9 

-23.0 

4.9 

9.2 

3 

1549.5 

+2.5 

-0.6 

21.5 

-18.1 

-19.6 

- 

-24.2 

4.6 

8.6 

4 

1553.2 

+3.5 

+0.9 

24.6 

-20.3 

-20.6 

-24.6 

- 

4.4 

8.2 

Table  1:  Packaged  MGF  (device  #7),  input  power=  - 9  dBm  (tunable  DBR  laser),  T~  20  0  C  with  currents 
of  26  mA  in  LA,  17  mA  in  SA  (input  side)  and  20  mA  for  grating  tuning.  The  transparency  current  of  each 
amplifier  for  each  channel  wavelength  is  also  given. 


unpolarized  transmission  spectrum  of  a  packaged  de¬ 
vice  is  shown  in  Figure  2.  The  channel  wavelengths 
are  spaced  by  3.7  nm  (see  Table  1).  Temperature  sen¬ 
sitivity  of  the  channel  positions  is  about  0.1  nm/K. 

A  tunable  DBR  laser  was  used  to  measure  the 
transparency  currents  of  the  long  (LA)  and  the  short 
(SA)  amplifier  at  each  channel  wavelength  (Table  1). 
Then  the  gain  and  gain  saturation  of  each  amplifier 
was  studied  on  a  packaged  device  by  measuring  the 
transmission  as  a  function  of  current  and  input  opti¬ 
cal  power  (with  the  other  amplifier  biased  at  trans¬ 
parency).  Results  are  shown  in  Figure  3  for  the  long 
amplifier  (LA)  with  channel  4.  The  solid  and  dashed 
curves  are  the  maximum  and  minimum  transmission 
(respectively)  when  the  input  polarization  is  varied. 
The  other  channels  and  SA  give  similar  results:  there 
seems  to  be  an  intrinsic  3  dB  loss  difference  between 
the  two  polarizations  and  gain  saturation  appears 
above  -5  dBm  input  optical  power. 

Fiber-to-fiber  losses  and  worst  case  crosstalk  (for 
any  polarization  combination)  for  a  packaged  MGF 
are  shown  in  Table  1.  Amplifier  currents  were  26 
mA  and  17  mA  (grating  tuning  current  =  20  mA), 
limited  by  the  onset  of  lasing,  and  the  input  power 
was  -9  dBm  (SA  side).  For  all  channels,  we  have  a 
fiber-to-fiber  gain  of  about  1  dB,  with  about  3  dB 
polarization  dependence.  The  crosstalk  is  between 
—  18  dB  and  —25  dB  in  all  cases.  We  estimated  that 
the  propagation  loss  on  the  chip  (including  the  butt- 
joints)  is  only  about  9  dB,  while  the  amplifiers  give 
about  16  dB  and  8  dB  gain  respectively.  This  means 
that  we  have  about  6.5  dB  coupling  loss  per  facet, 
which  is  high  and  may  be  due  in  part  to  the  use  of 
non-optimal  lenses. 

One  remaining  problem  is  the  origin  of  the  3  dB 
polarization  dependence  of  the  losses.  Part  of  this 
may  be  due  to  the  design  of  the  amplifier  sections, 
which  have  a  calculated  confinement  factor  of  0.296 
for  TE  and  0.231  for  TM.  And  even  in  the  gratings, 
free-carrier  loss  will  be  different  for  TE  and  TM  due 
to  different  field  distributions. 

In  conclusion,  we  have  designed,  fabricated  and 
packaged  a  four-channel  Multi-Grating  Filter  with 


polarization  independent  channel  wavelengths,  no 
fiber-to-fiber  loss  and  crosstalk  below  -18  dB. 

This  work  was  supported  in  part  by  the  RACE 
MWTN  (R2028)  project. 
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Abstract: 


Crosstalk  in  acousto-optic  switches  is  modulated  at  frequencies  related  to  that  of  the 
acoustic  waves  driving  the  device.  We  present  measurements  of  the  modulation,  a  brief 
discussion  of  its  causes  and  examples  of  its  effect  on  system  performance. 
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Introduction:  The  interaction  of  channels  in 

acoustooptic  switches  has  not  yet  been  fully 
characterized.  In  this  paper  we  report  observations  of 
modulated  crosstalk  in  these  devices,  and  show  that  two 
different  interacation  mechanisms,  both  associated  with 
the  traveling-wave  nature  of  the  acousto-optic 
interaction,  contribute  to  the  observed  modulation. 
This  crosstalk  can  critically  affect  the  performance  of 
the  switch  in  optical  wavelength-routing  crossconnects. 

The  acousto-optic  switch,  based  on 
wavelength-dependent  acousto-optic  polarization 
conversion,  is  being  investigated  as  a  wavelength 
routing  switch  in  WDM  networks,  primarily  because  it 
had  been  believed  to  offer  simultaneous  and 
independent  switching  of  many  closely-spaced 
wavelengths.  Recent  systems  tests  have  shown  that 
channel  independence  is  only  relative;  nearby  channels 
are,  on  the  average,  shifted  towards  each  other,  and 
switching  efficiency  is  decreased.1  We  have  shown 
recently  that  if  the  acousto-optic  polarization  converter 
is  redesigned  such  that  nearest  neighbor  channels  are 
switched  in  separate  sections  of  the  device,  these 
interactions  are  decreased,  substantially  restoring  the 
independence  of  the  channels.2 

A  related  interaction  between  channels,  which 
had  been  observed  in  the  past,  is  the  difference- 
frequency  modulation  of  the  signal  when  a  second, 
nearby  channels  is  selected.3  That  is,  when  two 
acoustic  frequencies,  and  £T2,  are  used 

simultaneously,  the  optical  output  is  found  to  be 
modulated  at  =  Qj  -  Q,2.  We  will  show  here  that 
modulation  occurs  at  the  acoustic  driving  frequencies, 
£2i,  and  at  2 £lly  as  well  as  at  die  sum  and  die  previously 
observed  difference  frequencies,  and  that  the 
dependence  of  modulation  depth  on  the  nature  of  the 
modulated  light  indicates  that  at  least  two  processes  are 
responsible.  Because  not  all  of  the  modulation  can  be 
eliminated,  this  effect  may  degrade  the  performance  of 
acousto-opdc  switches  when  large  numbers  of  channels 
are  selected  simultaneously. 

Sources  of  sum  and  difference  frequencies:  There 
are  at  least  two  sources  of  sum  and  difference  frequency 


modulation  in  die  AOTF,  beating  between  frequency 
shifted  components,  and  modulation  at  the  acoustic 
frequencies  of  die  polarization  conversion  efficiency. 

(a)  Frequency  shifts:  Beating  between  frequency- 
shifted  optical  components  was  commented  on  by 
Clioy,  et.  al.,  in  dieir  paper  describing  observation  of 
difference  frequencies.  When  two  acousuc  waves  are 
present,  each  produces  polarization  conversion  at  every 
wavelengdi.  Conversion  efficiency  is  determined  bodi 
by  die  acousdc  amplitude  and  by  the  extent  to  which  it 
reconciles  the  TE-TM  phase-mismatch.  When  phase 
matching  is  not  accurate,  conversion  efficiency  is  low; 
thus  die  device  is  highly  wavelengdi  selective. 
Additionally,  there  is  an  accompanying  opdcal 
frequency  shift4  equal  to  die  frequency  of  the  acoustic 
wave.  TE  and  TM  components  are  shifted  in  opposite 
directions,  and  the  direction  of  die  shift  depends  on  die 
reladve  direction  of  propagation  of  the  optical  and 
acousdc  waves.  Thus,  when  two  channels  are  switched 
simultaneously,  for  light  at  any  initial  opdcal 
wavelengdi  die  polarization-converted  light  has  two 
frequency-shifted  components,  one  for  each  of  the 
acousdc  waves;  these  beat  widi  each  other  to  produce 
sum  or  difference  frequencies,  depending  on  die 
relative  propagation  direction  of  the  two  acoustic 
waves.  We  expect  diat: 

•Such  beating  will  be  present  whedier  or  not  the  two 
acoustic  waves  overlap  in  die  device. 

•Beating  should  be  seen  only  to  die  extent  diat  die 
opdcal  source  is  coherent. 

•The  size  of  the  modulation  will  depend  on  the  reladve 
sizes  of  the  two  interfering  waves,  and  therefore  on 
bodi  die  conversion  efficiency  in  a  the  desired  channel 
and  diat  due  to  sidelobes  of  odier  channels 

In  addition,  if  polarization  isoladon  is  imperfect,  up 
and  down  conversion  leads  to  moduladon  at  die 
acousdc  driving  frequency  and  at  twice  that  frequency, 
since  die  polarization-converted  light  can  beat  with 
unconverted,  and  dierefore  unshifted  light  or  widi  light 
diat  has  been  shifted  in  die  opposite  direcdon.  Thus  we 
expect  moduladon  at  Q.  and  at  2f2,  even  when  only  one 
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acoustic  frequency  is  present,  and  at  sum  (or  difference) 
frequencies  where  only  difference  (or  sum)  would  be 
allowed  in  the  perfect  polarization  case. 

(b)  Modulation  of  conversion  efficiency:  A  second 
mechanism5'6  which  generates  sum  or  difference  a. 
frequencies  is  the  direct  modulation  of  polarization 
conversion  efficiency.  Although  the  acousto-optic 
interaction  leading  to  polarization  conversion  is  linear, 
the  dependence  of  the  polarization  conversion  on 
acoustic  amplitude  is  not.  Indeed,  the  AOTF  is  really  a 
directional  coupler,  a  notoriously  nonlinear  device, 
with  the  acoustic  wave  providing  the  coupling  between 
the  TE  and  TM  optical  waves,  and  with  a  response 
proportional  to  sin(icL),  where  k  is  the  coupling 
strength  and  L  the  length  of  the  device.  With  two 
acoustic  waves  present  in  the  device  simultaneously, 
the  coupling  coefficient  is  modulated  at  their  difference 
or  sum  frequency,  depending  of  relative  propagation 
direction  for  the  acoustic  waves,  leading  to  modulation 
of  the  polarization  conversion  efficiency  at  those 
frequencies.  When  the  polarization  conversion  is 
overdriven  (i.e.,  idL>Jt/2)  higher  harmonics  of  the 
difference  (or  sum)  will  also  be  generated.  We  can 
conclude  that' 

•Sum  and  difference  frequency  modulation  due  to 
nonlinearity  is  proportional  to  the  overlap  of  the 
acoustic  waves. 

•Modulation  does  not  require  a  coherent  light  source. 

The  polarization  conversion  efficiency  is  being 
modulated,  and  the  nature  of  the  optical  signal  is 
irrelevant 

•The  size  of  the  modulation  is  not  related  to  the 
sidelobe  intensity. 

Device  Configuration  and  Measurements:  We  have 
made  all  measurements  using  a  passband- flattened 
apodized*  AOTF,  with  an  acoustic  absorber  dividing 
the  interaction  region  in  half  and  two  IDTs,  which 
allow  each  half  of  the  device  to  be  addressed 
independently.  If  the  absorber  is  complete,  we  call  this 
an  apodized  device,  and  if  the  absorber  allows  some  of 
the  acoustic  wave  to  pass  through  it,  this  is  a  passband- 
flattened  AOTF.  Thus  we  can  look  at  copropagating 
acoustic  waves  with  total  overlap,  and 
counterpropagating  acoustic  waves  with  no  overlap  or 
with  small  overlap,  but  not  non-overlapping  acoustic 
waves  propagating  in  the  same  direction. 

In  order  to  characterize  the  sum  and  difference 
frequency  modulation,  input  light,  either  the  locally 
“white  light”  of  an  EDFA,  a  narrow  line  laser,  or  the 
laser  output  amplified  through  the  EDFA,  is  polarized, 
passes  through  the  AOTF,  and  then  through  a  second 
polarizer.  The  output  is  split,  with  half  going  to  an 


Fig.  I  Difference  frequency  modulation  of  bar  state 
output  (a)  at  500  kHz,  where  input  is  a  laser  modulated 
at  300  kHz  and  (b)  at  200  kHz,  with  unmodulated  input; 
with  wavelength  interaction  overdriven,  higher 
harmonics  are  produced.  Bar  state  average  rejection  is 
-16dB. _ _ _ 

optical  spectrum  analyzer  and  half  to  an  optical  RF 
spectrum  analyzer.  Thus  we  can  monitor  both  the 
optical  and  the  RF  spectrum  simultaneously.  In  an 
alternative  measurement,  part  of  the  optical  output  is 
detected  and  displayed  on  an  oscilloscope  (Fig.  1). 

Results:  For  all  measured  device  configurations  both 
sum  and  difference  frequency  modulations  were 
observed,  and  in  addition,  the  optical  output  was 
modulated  at  both  the  acoustic  driving  frequency,  £1 
and  at  2D.  Fig.  2  shows  a  typical  RF  spectrum.  In 
some,  but  not  all  cases  we  also  observed  modulation  at 
2(f2rfi2)  and  higher  harmonics.  Relative  modulation 
depths  depend  on  the  device  configuration,  the  nature 
of  the  optical  signal,  and  the  frequency  separation  of 
the  acoustic  waves.  The  data  indicate  that  both 
proposed  mechanisms  are  contributing.  With  a  laser 
input,  we  see  modulation.  Independent  of  acoustic 
wave  overlap,  as  predicted  from  beating  between 
frequency  shifted  components.  As  expected,  the  depth 
of  modulation  is  correlated  with  sidelobe  height  The 
modulation  depth  at  Q  and  2fl  could  be  reduced  almost 
to  zero  by  improving  polarization  isolation,  also  in 
agreement  with  the  above  arguments. 

However,  with  the  EDFA  alone  used  as  input  we 
expect  no  contribution  from  this  mechanism,  and 
indeed  the  observed  modulation  is  strongly  dependent 
on  both  acoustic  wave  overlap  and  relative  acoustic 
wave  propagation  direction.  Strong  modulation  is  seen 
only  for  the  case  where  both  acoustic  waves  propagate 
in  the  same  direction  in  the  same  section  of  the  device, 
and  only  the  difference  frequency  is  strong. 

The  worst-case  difference  frequency  modulation  is 
about  -10  dB.  Modulation  due  to  up  and  down 
conversion  can  be  minimized  by  increasing  peak 


power 


PD5-4 


200-1200  kHz 

174.5  -175.5  MHz 
349  -  351  MHz 


Fig.  2  RF  spectra  covering  three  frequency  ranges,  superimposed  to  show  simuftneous  generation  of  sum, 
difference,  acoustic  drive,  and  twice  acoustic  drive  frequencies.  Because  the  switch  is  “overdriven”  we  also  see 
twice  the  difference  frequency. 


conversion  efficiency,  flattening  the  passband,  and 
decreasing  sidelobes;  modulation  of  the  conversion 
efficiency  appears  to  be  intrinsic  to  the  device,  but  can 
be  limited  by  increasing  the  channel  spacing  to  channel 
width  ratio,  and  by  lengthening  the  device. 

Consequences:  The  presence  of  modulation  at  these 
frequencies  has  little  effect  on  the  signal  intended  to  be 
passed  through  the  filter.  Even  the  worst  modulation  of 
about  -10  dB  results  in  only  a  10%  variation  in  the 
intensity  of  the  signal.  However,  the  effect  on  crosstalk 
is  more  severe.  A  modulation  of  10%  of  the  signal 
implies  that  the  worse  case  crosstalk  will  be  at  least  -10 
dB.  While  switch  dilation  can  reduce  this  to  -20  dB, 
this  level  of  crosstalk  is  great  enough  to  limit  the 
usefulness  of  this  switch  for  certain  applications.  For 
example,  we  have  observed  a  power  penalty  for  high  bit 
rate  signals  when  four  close-spaced  channels  are 
switched  simultaneously.  Fig.  4  shows  BER 

BER 
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1.00E-03 
1.00E-04 
1.00E-05 
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i.ooe-08 
i.ooe-oe 
i.ooe-io 

-30.5  -29.5  -28.5  -27.5  -26.5  -25.5 

received  power 


Fig.  3  BER  measurements  made  at  2.48  Gb/s, 
comparing  behavior  of  acoustooptic  switch  with 
one  (1000),  three  (1 1 10),  or  four  (1111)  channels 
switched. 


measurements  made  at  2.48  Gb/s  with  1,3,  and  4 
channels  selected.  We  believe  that  the  increase  in  BER 
is  at  least  partly  attributable  to  time-dependent 
crosstalk. 
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Abstract 

Spatially  resolved  measurement  of  waveguide  birefringence  and  acoustic  interaction  strength  in 
Acousto-Optic  Tunable  Filters  is  demonstrated  using  short  acoustic  pulses. 


PD6-2 


Acoustic  Pulse  Measurements  of  Acousto-Optic  Tunable  Filter  Properties 
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Summary 

There  is  a  great  deal  of  interest  in  Acousto-Optic  Tunable  Filters  (AOTFs)  for  wavelength 
division  multiplexing  applications  because  of  their  fast  tuning  speeds,  wide  tuning  ranges  and  the 
ability  to  switch  multiple  wavelengths  simultaneously.'  One  practical  limitation  of  these  devices 
has  been  sidelobes  levels  which  are  often  much  larger  than  expected  by  theory.  This  is  the  case 
both  in  uniform  devices  as  well  as  those  where  apodization  techniques  to  reduce  sidelobes  are 
employed.  Trutna  el  al.2  explained  these  anomalous  sidelobes  as  being  the  result  of  longitudinal 
birefringence  variations  in  the  AOTF  waveguides.  These  variations  could  be  due  to  numerous 
sources  such  as  waveguide  width  and  doping  variations,  thermal  gradients,  substrate 
birefringence  variations  and  stress  induced  birefringence  variations.  While  it  was  useful  to  know 
the  general  source  of  high  sidelobes,  the  shape  of  the  filter  function  does  not  directly  yield  the 
location  of  the  problem,  or  help  with  finding  a  solution. 

Here  we  report  a  new  measurement  technique  which  gives  a  spatially  resolved  map  of  the 
birefringence  of  the  waveguide  as  well  as  the  strength  of  the  SAW  wave  and  the  acousto-optic 
interaction.  Additionally,  we  give  two  measurement  examples  which  demonstrate  the  usefulness 
of  the  technique  in  device  development.  The  basic  experimental  setup  is  shown  in  Figure  1. 
Polarized  light  is  launched  into  the  AOTF  and  the  output  passes  through  a  crossed  polarizer  and 
onto  a  detector.  The  acoustic  transducer  of  the  AOTF  is  driven  with  a  short  pulse  of  RF. 
Because  the  fractional  bandwidth  of  a  typical  AOTF  transducer  is  10%  or  more,  it  is  possible  to 
generate  acoustic  pulses  as  short  as  0. 1  ps.  Based  on  the  SAW  velocity  of  3750  m/s,  this 
corresponds  to  an  acoustic  pulse  length  of  0.375  mm  which  is  much  shorter  than  the  device.  We 
can  then  record  the  output  as  a  function  of  time  as  the  pulse  travels  down  the  device  length.  The 
resulting  output  is  a  spatial  map  of  the  acousto-optic  interaction  strength  along  the  length  of  the 
device  for  the  RF  frequency  used. 
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Fig  1 .  Basic  setup  for  pulsed  acoustic  measurements 
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Variations  in  the  acousto-optic  interaction  strength  can  be  due  to  two  causes: 

1)  Variations  in  the  acousto-optic  coupling  constant  k. 

2)  Variations  in  the  phase  mismatch  parameter,  Ap  (birefringence). 
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We  have  identified  two  measurement  modes  based  on  short  or  long  pulses  which  can  be  used  to 
sort  out  these  two  effects.  First,  let  us  examine  the  sensitivity  to  width  and  temperature 
variations  as  a  function  of  pulse  length.  Assuming  that  the  coupling  strength  and  phase 
mismatch  are  uniform  within  the  pulse,  the  converted  output  is  given  by: 


P  =- 


(k L)2 


(kL)2  + 


fAPL] 

2  J 


sin 


(kL)2  + 


fA|3L 


V 


(1) 


where  L  =  vax  is  the  pulse  length,  k  is  the  coupling  constant  and  A|3  is  the  phase  mismatch: 

2jtAn  2%fa 
X  v. 


Ap 


(2) 


For  the  typical  short  pulse  case  where  k L  « 7t/2,  the  FWHM  of  equation  (1)  in  terms  of  Af3  is: 

aD  5.6  5.6 

^ FWHM  ~  i  ~  (KL  «  1)  (3) 

JLm  V„T 


Using  measured  dependencies  of  brifringence  on  waveguide  width  and  temperature,  the  required 
waveguide  width  and  temperature  variations  to  reduce  the  output  by  half  are  then: 


8w,/2  - 


1.15pm 

-t(iu) 


STm~ 


5.35°  C 

x(pj) 


(4) 


If  a  very  short  pulse  is  used  (say  0. 1  ps),  the  output  is  found  to  be  very  insensitive  to  realistic 
birefringence  variations.  For  example,  the  waveguide  width  would  have  to  vary  by  11.5  pm  or 
the  temperature  by  53 °C  to  reduce  the  output  by  half.  Given  that  Ap  is  effectively  0  in  the  short 
pulse  regime,  we  find  that  the  output  becomes: 

Poul (short pulse)  =  sin2(tcL) «  (kL)2  ~  Pa  (5) 


Thus  the  output  becomes  a  map  of  the  square  of  the  coupling  constant  only  which,  if  we  assume 
a  uniform  acousto-optic  coefficient,  is  proportional  to  the  acoustic  power.  An  example  of  this 
measurement  is  shown  in  Figure  2.  This  data  shows  that  there  is  not  undo  attenuation  or  walk- 
off  of  the  acoustic  wave. 


Equivalent  Distance  (mm) 


Fig  2.  Short  pulse  measurement  results  showing  acoustic  power  uniformity 


In  the  long  pulse  regime  (e.g.  2  ps)  there  is  reduced  spatial  resolution,  but  the  output  is  sensitive 
to  phase  mismatch  as  well  as  the  coupling  constant.  By  taking  a  series  of  measurements  at 
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different  acoustic  frequencies,  we  can  build  up  a  contour  map  as  shown  in  Figure  3.  In  this  plot, 
we  have  made  time  (which  is  equivalent  to  position  along  the  waveguide)  the  x  axis.  The  y-axis 
is  the  acoustic  frequency.  In  an  ideal  device  with  an  on  axis  transducer,  the  plot  should  show  a 
uniform  horizontal  ridge  centered  on  the  phasematching  frequency.  Birefringence  variations  will 
be  seen  as  curvature  of  that  ridge,  and  acoustic  power  variations  will  be  seen  as  variations  of  the 
ridge  height.  The  device  measured  in  Figure3  is  shown  on  the  right.  It  is  a  standard  single 
section  AOTF  with  the  transducer  offset  laterally  and  rotated  such  that  the  acoustic  beamwalks 
on  and  then  off  the  waveguide  with  the  aim  of  providing  a  simple  form  of  amplitude  apodization. 
In  the  resulting  contour  plot,  we  can  see  both  the  amplitude  apodization  as  well  as  a  small 
amount  of  curvature  indicating  a  birefringence  nonuniformity.  The  righthand  axis  of  the  plot 
calibrated  the  birefringence  offset  in  terms  of  equivalent  waveguide  width  variations. 

We  have  already  used  this  technique  to  solve  a  packaging  problem  where  stress  induced 
birefringence  varied  with  time.  In  that  case,  the  shape  of  the  filter  function  varied  as  the 
temperature  of  the  package  was  cycled,  with  a  very  long  recovery  time.  Using  these  techniques, 
it  was  possible  to  identify  where  the  chip  was  epoxied  to  the  package  and  where  it  was  free  to  be 
stretched.  The  ability  to  measure  the  acoustic  power  as  a  function  of  distance  should  also  be 
useful  in  designing  amplitude  apodization  schemes.  Previously,  the  only  was  to  evaluate  these 
ideas  was  to  use  a  SAW  probe,  which  requires  a  highly  specialized  measurement  setup. 

In  summary,  we  have  described  a  new  measurement  technique  for  AOTFs  which  provides 
spatially  resolved  data  on  the  birefringence  and  acoustic  power  non-uniformities  in  the 
waveguide.  This  technique  is  simple  to  implement  and  requires  no  special  equipment  ither  than 
a  pulsed  RF  source.  In  addition  to  various  AOTF  applications,  this  tehcnique  may  be  valuable 
for  evaluating  waveguide  or  package  stress  problems  in  other  LiNb03  integrated  optical  devices. 
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Fig  3.  Long  pulse  measurements  at  various  frequencies  showing  waveguide  uniformity 
and  amplitude  apodization  in  and  AOTF  with  an  offset  rotated  transducer  (right). 
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Abstract 

Basic  design  rule  for  highly-efficient  operation  of  1.3-/im  InP-based  strained  MQW  lasers 
at  elevated  temperatures  is  derived  on  the  basis  of  a  self-consistent  numerical  approach. 
The  validity  of  the  design  rule  is  experimentally  demonstrated  by  the  excellent  perfor¬ 
mance  (77=0.55  W/A  and  Pmax  >  35  mW  at  363  K)  of  fabricated  lasers. 
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Basic  Design  Rule  of  1.3-fim  InP-Based  Strained  MQW  Lasers 
for  High  Temperature  Operation 

S.  Seki,  H.  Oohashi,  H.  Sugiura,  T.  Hirono,  and  K.  Yokoyama 
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InP-based  strained-layer  (SL)  multiple-quantum-well  (MQW)  lasers  operating  at  a 
wavelength  of  1.3  ym  are  of  much  interest  due  to  their  applications  to  fiber-in-the-loop 
(FITL)  and  fiber-to-the-home  (FTTH)  systems  [1],  [2].  For  these  applications,  high 
efficiency  and  high  output  power  are  required  as  well  as  low  threshold  current  over  a 
wide  range  of  operation  temperatures.  In  this  paper,  we  derive  a  basic  design  rule  for 
highly-efficient  operation  of  InP-based  SL-MQW  lasers  at  elevated  temperatures  on  the 
basis  of  a  self-consistent  numerical  approach. 

The  band  offset  in  InP-based  SL-MQWs  is  significantly  changed  under  high- 
injection  conditions  [3],  [4],  To  take  into  account  the  changes  in  band  offset  due 
to  electrostatic  band-profile  deformation,  the  band  structures  were  obtained  by  self- 
consistently  solving  the  Poisson  equation,  the  scalar  effective-mass  equation  for  the 
conduction  band,  and  the  multi-band  effective-mass  equation  for  the  valence  band.  To 
properly  include  the  strain-dependent  coupling,  a  6  x  6  Luttinger-Kohn  Hamiltonian 
was  used  for  the  valence-band  Hamiltonian.  The  conduction-band  Hamiltonian  was 
assumed  to  be  parabolic  in  k-space. 

In  this  study,  we  focus  on  two  design  parameters:  (i)  the  ratio  of  radiative  current 
component,  47  / 47*  >  and  (ii)  the  threshold  carrier  density.  Here,  J^^is  the  radia¬ 
tive  current  component  resulting  from  the  transition  between  the  ground  states  in  the 
conduction  and  valence  bands,  and  4aal  is  the  total  radiative  current  density.  Hence, 
4-d74dal  is  a  good  measure  of  carrier  confinement  in  the  QWs. 

We  analyzed  an  MQW  structure  with  compressively-strained  InAsyPi_y  (y  =  0.52) 
[5]  as  the  active  region.  The  amount  of  strain  and  the  QW  thickness  were  adjusted  to 
be  1.5%  and  6  nm  [6].  The  MQW  stack,  which  incorporates  three  QWs  with  barriers 
of  10-nm  InGaAsP,  is  sandwiched  between  50-nm  InGaAsP  separate  confinement  het¬ 
erostructure  (SCH)  layers.  To  clarify  the  influence  of  carrier  spillover,  three  different 
compositions  of  InGaAsP  were  used  for  the  SCH  layers:  Ag  =  1.1,  1.05  and  1.0  ym  (re¬ 
ferred  to  hereafter  as  Ql.l,  Ql.05,  and  Q1.0,  respectively).  The  composition  of  barrier 
layers  is  the  same  as  that  of  SCH  material  in  each  case. 

The  physical  meaning  of  47  / 4aal  *s  more  clearly  indicated  in  Fig.  1,  which  shows 
the  spontaneous  emission  spectra  for  SL-MQWs  with  Ql.l  barrier  and  SCH  layers  when 
the  mode  gain  T  Gp  is  30  cm-1.  Here,  T  is  the  optical  confinement  factor  and  Gp  is 
the  peak  value  of  optical  gain.  The  shaded  area  corresponds  to  the  radiative  current 
component  4a1-  It  is  noted  that  the  peak  corresponding  to  the  band-gap  energy  of  the 
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barrier  and  SCH  material  (i.e.,  about  1.1  eV)  is  notably  enhanced,  which  agrees  well 
with  experimentally  observed  spectra  changes  [7].  These  results  clearly  indicate  that 
it  is  necessary  to  take  into  account  electrostatic  band-profile  deformation  to  accurately 
determine  the  design  parameter  J^1/ at  elevated  temperatures. 

Figure  2  shows  the  variations  in  / *^radal  relative  to  mode  gain,  using  the  com¬ 

position  of  the  barrier  and  SCH  material  as  a  parameter.  The  results  for  two  different 
temperatures  are  compared.  Note  that  the  band-gap  energy  of  the  barrier  and  SCH  ma¬ 
terial  has  a  pronounced  influence  on  4a1  / ^radal  elevated  temperatures.  These  results 
suggest  that  the  efficiency  of  lasers  at  elevated  temperatures  is  significantly  enhanced  by 
using  the  barrier  and  SCH  materials  with  larger  band-gap  energy.  However,  increasing 
band-gap  energy  of  these  materials  could  significantly  reduce  the  optical  confinement 
factor,  leading  to  an  increase  in  threshold  carrier  density  (iVth).  Hence,  it  is  necessary 
to  make  a  trade-off  between  4‘d1  /jtota1  ancj  jyth>  The  variations  in  Nt h  with  changes 
in  mode  gain  are  plotted  in  Fig.  3,  using  the  composition  of  barrier  and  SCH  layers  as 
a  parameter.  Over  a  wide  range  of  mode  gain  (3-10  cm-1 /well),  Ql.05  gives  the  lowest 
Nth.  regardless  of  temperature,  while  Ql.O  increases  Nt h-  These  results  indicate  that, 
to  maximize  device  performance  for  a  required  mode  gain  at  elevated  temperatures, 
SL-MQW  structures  should  be  designed  so  as  to  maximize  J^1/ Jrl°^al  without  paying 
much  penalty  in  Nth  ■ 

Applying  the  derived  design  rule,  we  fabricated  SL-MQW  Fabry-Perot  buried 
heterostructure  (BH)  lasers.  The  SL-MQW  stack,  which  consists  of  six  InAs^P^j, 
(y  =  0.52)  QWs  with  1.5%  compression,  was  grown  by  chemical  beam  epitaxy  [8],  In 
accordance  with  the  results  shown  in  Figs.  2  and  3,  we  chose  Ql.05  for  the  barrier  and 
SCH  layers.  Figure  4  shows  the  CW  light-current  characteristics  of  SL-MQW  BH  lasers. 
The  laser  cavity  is  300-/xm  long  with  a  high  reflection  coating  (Rr  =  96%)  on  the  rear 
facet.  The  threshold  current  is  6.5  mA  at  303  K  and  the  characteristic  temperature 
(T0)  of  threshold  current  is  56  K  in  the  temperature  range  of  303-373  K.  Note  that  very 
high  efficiency  of  0.55  W/A  is  obtained  at  363  K.  This  makes  it  possible  to  reduce  the 
drive  current  required  for  high  power  operation  at  elevated  temperatures.  At  363  K, 
the  output  power  is  30  mW  at  drive  currents  as  low  as  100  mA.  This  output  power  at 
363  K  is  more  than  one  and  half  times  as  large  as  that  when  Ql.l  or  Ql.O  was  used  for 
the  barrier  and  SCH  layers  under  the  same  current  condition.  The  maximum  possible 
power  at  363  K  exceeds  35  mW  by  further  increase  in  drive  current.  These  results 
clearly  indicate  the  validity  of  the  derived  design  rule. 

In  summary,  we  have  derived  a  basic  design  rule  of  1.3-^m  InP-based  SL-MQW 
lasers  for  high  temperature  operation.  The  design  rule,  which  was  obtained  on  the 
basis  of  a  self-consistent  numerical  approach  including  the  Poisson  equation  and  the 
effective-mass  equations  for  the  conduction  and  valence  bands,  is  to  maximize  the  ratio 
°f  ^rad1  / ^rad*-1  f°r  a  required  mode  gain  without  paying  much  penalty  in  iVth.  Following 
this  rule,  high-efficiency  (0.55  W/A  at  363  K)  and  high-power  (over  35  mW  at  363  K) 
InP-based  SL-MQW  BH  lasers  have  been  fabricated. 
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Fig.  1  Calculated  spontaneous  emission  Fig.  2  Current  component  ratio  plotted  as  a 
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MQW  lasers.  of  the  barrier  and  SCH  material  as  a 

parameter. 


Current  (mA) 

Fig.  3  Threshold  carrier  density  plotted  as  a  Fig.  4  CW  light-current  characteristics  of  a 
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Abstract 

We  outline  and  examine  numerically  an  efficient  full  vectorial  beam  propagation 
algorithm  based  on  the  alternating  direction  implicit  (ADI)  finite  difference 
method. 
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Introduction:  While  relaxation  procedures  for  determining  the  vector  modes  of  waveguides  can  be 
immediately  converted  into  imaginary  distance  propagation  methods, [1]  full  vector  beam  propagation 
techniques  in  real  space  are  intrinsically  unstable[2],  Further  as  alternating  directional  methods  can  only  be 
applied  to  a  particularly  divergent  form  of  the  vector  wave  equation  (although  they  are  easily  adapted  to  the 
quasivectorial  approximation)^],  the  computational  effort  required  in  current  propagation  programs  is  often 
excessive.  Direct  finite  difference  methods,  while  rapid,  are  only  conditionally  stable,  and  thus  may  require 
small  step  sizes  for  certain  refractive  index  profiles[4]. 

In  this  paper,  we  consider  a  novel  alternating  directional  implicit  (ADI)  procedure,  first  proposed  in  [5],  for 
the  full  vector  wave  equation.  This  method  is  shown  to  be  rapid  and  simple  to  implement  numerically.  We 
then  apply  the  procedure  to  the  test  case  of  a  longitudinally  varying  semiconductor  rib  waveguide  Y- 
junction,  obtaining  results  which  yield  physical  insight  regarding  the  mutual  coupling  between  the  two 
polarizations. 

Numerical  Method:  We  consider  the  vector  wave  equation,  which  can  be  written  in  the  compact  form 


322 


e  =  *2„2e  +  *  (*2  —  „2)  £  + 
a  00a  o  '  0'  a 


with  a,  p  €  {*,  y}  . 


alogn2  3  a2logn2 

3p  8a  8a8p 

H  can  be  immediately  separated  into  a  term  containing  only  x  derivative  and  a  second  term  containing  y 
derivative.  While  such  a  procedure  immediately  enables  the  construction  of  a  ADI  method,  the  presence  of 
second-order  cross-derivatives  of  the  refractive  index  gives  rise  to  a  high  degree  of  instability  in  the  resulting 
propagation  algorithms.  Accordingly,  we  instead  seek  a  ADI  procedure  for  the  unmodified  Eq.(1). 


If  the  off-diagonal  terms  in  the  operator  H  are  expressed  as 
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Introducing  the  paraxial  approximation  e±  (x,  y,  z)  =  ^(x,  y,  z)  <•  Vo1  the  formal  solution  to  Eq.(1 )  can 
be  written  as 
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While  the  two  leftmost  exponential  operators  may  be  evaluated  using  standard  splitting  procedures,  a 
coordinate  transformation  is  required  to  apply  the  ADI  method  to  the  third  operator[6].  After  discretization 
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on  a  two-dimensional  uniform  grid  with  Nx  x  Ny  points  yields 

M  >P'5  1  =  — —  p,  °2  “3  +  1,  j  +  1  y  +  1,  j  ~  1  yi  ~  1,  /  +  1 
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with  a1=(n2i+1  j+1  /n2j+1j)-1,  a2=(-n2i+1  j_1  /n2i+1  j)+1 ,  a3=(-n2j_i  j+1  /n2^  j)+1 ,  and  a4=(n2i.1j.1  /n2^  j)-1  and 
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in  which  b1=(n2i+1j+1  /n2iij+1)-1,  b2=(-n2M  j+1  /n2i  j+i)+1,  b3=(-n2i+1j.1  /n2i(j.i)+1,  and  b4=(n2M  H  /n2ij.1)-1. 
The  superscript  T  denotes  transpose. 

In  practice,  we  first  apply  a  one-dimensional  finite  difference  method  to  each  of  the  downward  sloping  and 
then  to  the  upward  sloping  diagonals  of  the  two-dimensional  grid  of  electric  field  values  shown  in  the 
following  schematic  diagram. 


Fig.  1  Field  Values  Coupled  by  the  ADI  Method 

As  the  terms  in  M^y  and  My^  along  the  upward  and  downward  diagonals  are  decoupled,  the  ADI  may 
be  employed.  For  the  Crank-Nicholson  finite  difference  method,  this  results  in  the  solution  of  2(NX+Ny) 
pentadiagonal  equation  systems.  Such  a  procedure  is  rapid  and  relatively  simple  to  implement  in 
comparison  with  sparse  matrix  inversion  techniques. 

Results:  To  illustrate  the  computational  method  of  the  proceeding  section,  we  consider  a  Y-junction  formed 
from  two  rib  waveguides  of  height  1.1  pm  and  width  W=2  pm  etched  into  a  1.3  pm  thick  epilayer  with 
refractive  index  n1=3.44.  The  cladding  layer  has  a  refractive  index  of  n^.34  while  the  light  wavelength  is 
set  to  the  \=1 .55  pm.  The  reference  refractive  index  n0  is  set  here  to  the  average  of  and  n2.  The  refractive 
index  profile  of  the  Y  junction  is  the  union  of  the  two  waveguide  profiles  centered  at  ±x(z)  where  x(z)=2[1  - 
cos(7rz/1 000pm)].  The  computational  window  dimensions  are  given  by  Lx  *  Ly=  16  *  8  pm2,  the  number  of 
transverse  grid  points  is  Nx  *  Ny  =64  *  64,  the  propagation  step  length  is  0.1pm  and  a  tapered  absorbing 
region  is  present  over  the  8  grid  points  closest  to  the  computational  window  boundaries. 

In  our  full  vectorial  beam  propagation  analysis,  the  waveguide  is  excited  with  the  lowest-order  TE  mode  (VFX, 
'Fy)  at  the  Y-junction  input.  The  intensities  of  both  ¥x  and  SFy  along  the  cross-section  located  at  the  at  the 
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etching  stop  layer  is  displayed  as  a  function  of  propagation  distance  in  Figs.  2  and  3,  respectively. 


Fig.2  Intensity  of  along  the  etching  stop  layer 
as  a  function  of  propagation  distance 


Fig.3  Intensity  ofH'y  along  the  etching  stop  layer 
as  a  function  of  propagation  distance 


We  observe  from  the  above  figures  that  the  divergence  associated  with  the  high-order  modes  is  suppressed 
by  the  absorber,  as  predicted  in  [5].  Note  the  periodic  transfer  of  electric  field  intensity  between  the  TE  and 
TM  field  components.  Calculations  performed  with  the  transparent  boundary  conditions  display  additional 
boundary  reflections  in  comparison  to  the  results  of  the  figures. 

Conclusion?:  We  have  demonstrated  a  rapid  yet  accurate  alternating  directional  procedure  for  the  vector 
wave  equation  which  diverges  relatively  slowly  compared  to  previous  formalisms.  Our  rib  waveguide  Y- 
junction  test  case  indicates  that  the  method  is  easily  stabilized  by  a  properly  tailored  absorber  and  thus  can 
be  employed  to  study  polarization  evolution  in  a  wide  variety  of  components. 
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1.  Introduction 

Ball  lenses  are  important  for  many  applications.  For  example,  ball  lenses  can  be 
used  to  match  the  mode  of  a  laser  diode  (LD)  to  a  single  mode  fiber  (SMF),  essential  for 
low-loss,  high  bit  rate  communication  systems.  Modeling  the  propagation  of  LD  light 
through  a  ball  lens  presents  a  challenge  due  to  the  large  angular  divergence  of  the  LD 
field  (typically  >20°  HWHM)  and  the  subsequent  significant  effect  of  spherical 
aberration.  Accurately  calculating  the  reflected  power  is  also  difficult,  but  essential, 
since  reflections  as  small  as  -30  dB  can  destabilize  the  LD.  [1] 

A  full-wave  analysis  of  this  system  using,  e.g.,  a  finite-difference  time-domain 
method  is  not  practical  because  of  the  size  of  the  ball  lens,  typically  hundreds  of 
wavelengths  in  diameter.  Approximate  scalar  methods  [2]  can  give  good  results  in  some 
cases,  but  fail  to  calculate  reflected  power  and  miss  polarization  effects  entirely. 

Our  approach  exploits  the  fact  that  the  scattering  of  an  arbitrary  electromagnetic 
beam  from  a  sphere  is  an  exactly  solvable  problem.  The  scattering  of  a  plane  wave  from 
a  sphere  is  a  classical  problem  which  was  solved  by  Mie  in  1908  [3].  More  recently, 
various  workers  have  considered  the  scattering  of  a  Gaussian  beam  from  a  sphere  [4] 
and  its  numerical  implementation  for  other  applications.  To  our  knowledge,  this  is  the 
first  time  this  approach  has  been  applied  to  a  problem  in  optical  design.  We  are  able  to 
calculate  reflection  and  transmission  accurately  with  modest  computational  effort. 

2.  Method 

We  assume  the  ball  lens  is  a  perfect,  homogeneous  sphere,  and  that  departures 
from  ideality  can  be  modeled  by  a  suitable  average  over  lens  radius  or  wavelength.  An 
LD  or  SMF  mode  is  modeled  as  an  approximately  linearly  (x)  polarized  (possibly 
astigmatic)  Gaussian  beam  of  vacuum  wavelength  l .  In  the  waist  plane  z  =  -z0  it  has 
the  form 

Ex{x,y,z  =  -z0 )  =  E0  exp(- (x  -  x0  )2  /  a]  -  (y  -  y0 )2/o] ) , 

where  xq,  yo  are  offsets  from  the  optical  (z)  axis.  We  need  to  know  the  field  everywhere 
outside  the  lens,  which  can  be  written  as  the  sum  of  the  incident  field  and  a  scattered 
field: 
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The  scattering  problem  is  solved  by  a  vector  multipole  expansion  of  the  incident  (LD  or 
SMF),  scattered  and  internal  fields.  Details  of  this  expansion,  which  involve  some  subtle 
convergence  issues  for  wavelength  scale  LD  spot  sizes,  will  be  presented  at  the 
conference.  Satisfying  the  boundary  conditions  at  the  surface  of  the  lens  (the  continuity 
of  r  x  E  and  f  x  H)  determines  the  field  everywhere.  In  practice  we  employ  a  T-matrix 
method  to  calculate  the  expansion  coefficients  of  the  field  [3].  We  require  typically 
-1000  terms  in  the  expansion  for  a  convergent  result.  Once  the  field  is  known  we  can 
calculate  overlap  integrals  of  interest.  For  example,  the  power  reflected  back  to  an 

optical  fiber  is  R  =  jj  E*„c  •  EJCa,  dxdy,  where  the  integral  is  over  the  plane  of  incidence 
z  =  -z0. 

3.  Results 

To  validate  our  method,  we  have  calculated  the  back  reflection  to  a  SMF  of  mode 
size  4.5  pm  (half-width  at  1  /e2  intensity)  at  various  lens-fiber  distances  for  two  different 
ball  lenses:  250  pm  sapphire  (n  =  1.75)  and  300  pm  BK-7  (n  =  1.5).  In  all  cases  the  fiber  is 

on  the  optical  axis.  The  light  wavelength  is  1.3  pm.  The  back  reflections  were  also 
measured  experimentally.  These  results  are  shown  in  Fig.  1;  it  is  apparent  that  the 
calculation  agrees  with  experiment  in  most  cases  to  better  than  1  dB. 
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Fig.  1.  (a)  Back  reflection  to  SMF  from  (a)  250  mm  radius  sapphire  lens;  (b)  300  mm 
radius  BK-7  lens 


The  SMF  results  give  us  confidence  to  calculate  the  back  reflection  to  an  LD, 
where  measurements  are  much  more  difficult  to  make.  In  Fig.  2  we  show  the  back 

reflection  to  an  LD  of  symmetric  spot  size  (0.85  pm  x  1.22  pm)2  as  the  LD-lens  distance 
is  varied.  For  comparison,  we  also  show  the  result  obtained  from  a  symmetric 

(1.22  pm)2  spot.  Note  the  nonmonotonic  behavior,  probably  due  to  high-order 
resonances  in  the  spherical  lens. 
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Fig.  2  Back  reflected  power  to  a  LD  for  a  300  mm  radius  ball  lens. 


Our  results  clearly  indicate  that  we  are  able  to  predict  accurately  the  reflection 
from  a  ball  lens  to  a  SMF.  We  believe  our  method  can  accurately  predict  the  back 
reflection  for  the  LD,  which  is  very  difficult  to  make  reliably  using  any  other  existing 
method.  This  technique  will  additionally  allow  us  to  study  the  effect  of  misalignments, 
and  the  extension  of  the  method  to  antireflection  coated  lenses  is  straightforward. 
Calculations  of  coupling  to  optical  fibers  will  also  be  forthcoming. 
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summary 

We  calculate  reflection  from  ball  lenses  by  exactly  solving  Maxwell's  equations  for  the 
scattering  of  a  beam  from  a  dielectric  sphere.  Our  results  are  validated  by  experiment. 
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